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LOCAL CALIBRATIONS FOR MINIMIZERS OF THE MUMFORD-SHAH 
FUNCTIONAL WITH A TRIPLE JUNCTION 

Maria Giovanna Mora 

Abstract. We prove that, if u is a function satisfying all Euler conditions for the 
Mumford-Shah functional and the discontinuity set of u is given by three line segments 
meeting at the origin with equal angles, then there exists a neighbourhood U of the 
origin such that u is a minimizer of the Mumford-Shah functional on U with respect 
to its own boundary conditions on dU . The proof is obtained by using the calibration 
method. 



1 Introduction 

The Mumford-Shah functional was proposed in ll^] to approach image segmentation problems and it can 
be written, in the "homogeneous" version in dimension two, as 

iVuix^^dxdy + n'iSu), (1.1) 

where Q is a bounded open subset of R 2 with a Lipschitz boundary, Ji 1 is the one-dimensional Hausdorff 
measure, u is the unknown function in the space SBV(fl) of special functions of bounded variation in 
f2, S u is the set of essential discontinuity points of u, while Vu denotes its approximate gradient (see 
1). 

This paper deals with local minimizers of (1.1) with given boundary values. More precisely, we say 
that u is a Dirichlet minimizer of (1.1) in Q if u belongs to SBV(fl) and satisfies the inequality 

/ \Vu(x,y)\ 2 dxdy + n\S u ) < [ \Vv{x,y)\ 2 dxdy + H 1 {S V ) 
Jn Jn 

for every v £ SBV(ft) with the same trace as u on dfl. 

Con sidering different classes of infinitesimal variations, one can show that, if u is a Dirichlet minimizer 
of ( |l.l|) in f2 , then the following equilibrium conditions (which can be globally called the Euler conditions 
for (Ml)) are satisfied (see [| and @): 

• u is harmonic on O \ S u ; 

• the normal derivative of u vanishes on both sides of S u , where S u is a regular curve; 
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• the curvature of S u (where defined) is equal to the difference of the squares of the tangential 
derivatives of u on both sides of S u ; 

• if S u is locally the union of finitely many regular arcs, then S u can present only two kinds of 
singularities: either a regular arc ending at some point, the so-called "crack-tip", or three regular 
arcs meeting with equal angles of 2w/3, the so-called "triple junction". 

However, since the functional (P-.lj) is not convex, the Euler conditions are not sufficient for the Dirichlet 
minimality of u. 

In |H it has been proved that, if S u is an analytic curve connecting two points of dfl (hence, S u 
has no singular points), then the Euler conditions are also sufficient for the Dirichlet minimality in small 
domains. In this paper we prove that, if S u is given by three line segments meeting at the origin with 
equal angles of 2tt/S (i.e., S u is a rectilinear triple junction), the same conclusion holds; in other words, 
for every (x, y) £ £1, there is an open neighbourhood U of (x, y) such that u is a Dirichlet minimizer of 



(1.1) in U. Since for (x,y) ^ (0,0) this fact follows from the result in J9|, the interesting case is when 
we restrict the functional to a neighbourhood of the triple point (0, 0) . 
The precise statement of the result is the following. 

Theorem 1.1 Let SI := _B(0, 1) be the open disc in R 2 with radius 1 centred at the origin, and let 
(Ao, Ai, A%) be the partition of defined as follows: 

2 2 1 
(r cos 9,r sin 9) G SI : < r < 1, -tt(2 - i) < 6 < -tt(3 - i) \ Vi = 0,l,2. 

3 3 1 

Let Sij := Ai D Aj for every i < j . Let Ui G C 2 {Ai) be a harmonic function in A4, satisfying the 
Neumann conditions on dAi n SI and such that |Vu»| = |Vuj| on Sij for every i < j . If u is the 
function in SBV{ST) defined by u := Uj a.e. in each Ai and t*o(0, 0) < ui(0, 0) < 1*2(0,0), then there 
exists a neighbourhood U of the origin such that u is a Dirichlet minimizer in U of the Mumford-Shah 
functional. 

This theorem generalizes the result of Example 4 in jlj , where the functions 11% were three distinct 
constants. The proof is obtained by the calibration method adapted in (lj to the functional (1.1). We 



construct an explicit calibration for u in a cylinder C/xR, where U is a suitable neighbourhood of 
(0,0). The symmetry due to the 27r/3-angles is exploited in the whole construction of the calibration; 
in particular, it allows to deduce from the other Euler conditions that each m must be either symmetric 
or antisymmetric with respect to the bisecting line of Ai and then, it can be harmonically extended to 
a neighbourhood of the origin, cut by a half-line in the antisymmetric case. Around the graph of it , the 
calibration is obtained using the gradient field of a family of harmonic functions, whose graphs fibrate a 
neighbourhood of the graph of u; this technique reminds the classical method of the Weierstrass fields, 
where the minimality of a candidate u is proved by constructing a suitable slope field, starting from a 
family of solutions of the Euler equation, whose graphs fibrate a neighbourhood of the graph of u . 

The assumption of C 2 -regularity for does not seem too restrictive: indeed, by the regularity results 
for elliptic problems in non-smooth domains (see [Q), it follows that Uj belongs at least to C 1 (Ai), since 
Ui solves the Laplace equation with Neumann boundary conditions on a sector of angle 27r/3. Moreover, 
since Ui is either symmetric or antisymmetric with respect to the bisecting line of Ai , one can see Ui as 
a solution of the Laplace equation on a 7r/3-sector with Neumann boundary conditions or respectively 
mixed boundary conditions. By the regularity results in H, it turns out that in the first case Uj belongs to 
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C 2 (Ai) , while in the second one Ui can be written as Ui(r, 9) — Ui(r, 9) + cr 3 / 2 cos |6* , with Ui £ C 2 (Ai) 
and c£t. So, only the function r 3 / 2 cos |6> is not recovered by our theorem. 

The case where S u is given by three regular curves (not necessarily rectilinear) meeting at a point 
with 27r/3-angles, is at the moment an open problem and it does not seem to be achievable with a plain 
arrangement of the calibration used for the rectilinear case, essentially because of the lack of symmetry 
properties. 

The paper is organized as follows. I n Se ction 2 we recall the main result of Q, while Sections 3 



7 are devoted to the proof of Theorem 1.1: in Section 3 we construct a calibration tp in the case Ui 
symmetric and we prove that ip satisfies conditions (a), (b), (c), and (e) (see the definition of calibration 
in Section 2); in Sections 4 and 5 we show some estimates, which will be useful in Section 6 to prove 
condition (d); finally, in Section 7 we adapt the calibration to the antisymmetric case. 

Acknowledgements. The author wish to thank Gianni Dal Maso for many interesting discussions on 
the subject of this work and for some helpful suggestions about the writing of this paper. 



2 Preliminary results 

Let f2 be an open subset of R 2 with a Lipschitz boundary. If u is a function in SBV(Vl) , for every 
{xqiDo) G fl one can define 

u + (xo,yo) := aplimsup u(x,y), u~(xo,yo):= apliminf u(x,y). 

We recall that S u — {(x,y) £ ft : u~(x,y) < u + (x,y)} and that for Ti}-a,.e. (xo,yo) £ S u there exists a 
(uniquely defined) unit vector v u {xo,yo) (which is normal to S u in an approximate sense) such that 



lim 



/ \u(x,y) - ^(xo^o^dxdy = 0, 

JB±(xa,Vn) 



r^0+ £ 2 (B±(x Q ,y )) Jb±{ X0 , V0 ) 

where B^(xQ,yo) is the intersection of the ball of radius r centred at (xo,j/o) with the half- plane 
{(x, y) £ R 2 : ±(x — xq, y — yo) ■ v u (xq, yo) > 0} . For more details see j|. 

For every vector field <p : JlxR^R 2 xR we define the maps ip xy : fixR R 2 and ip 2 : Hxl-tl 

by 

tp(x, y, z) = (ip xy (x, y, z), <p z (x, y, z)). 

We shall consider the collection T of all piecewise C 1 vector fields (/j:fixR^R 2 xR with the following 
property: there exist a finite family (£7j)iei of pairwise disjoint open subsets of f2xR with Lipschitz 
boundary whose closures cover f2xR, and a family (ipi)i^i of vector fields in C 1 (J7i,R 2 xR) such that tp 
agrees at any point with one of the tpi . 

Let u £ SBV(Q) . A calibration for u is a bounded vector field tp £ T satisfying the following 
properties: 

(a) div^ = in the sense of distributions in f2xR; 

(b) \<p xy (x, y, z)\ 2 < 4(p z (x,y,z) at every continuity point {x,y,z) of ip; 

(c) tp xy (x,y,u(x,y)) = 2Vu(x,y) and ip z (x, y, u(x, y)) = \Wu(x, y)\ 2 for a.e. (x, y) £ SI \ S u ; 
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(d) 



r*2 



(p xy (x,y,z)dz 



Hi 

r u + (x,y) 



< 1 for every {x,y) £ Q, and for every t\,t 2 G R; 



(e) / ip xy (x,y,z) dz = i> u (x 7 y) tor T(}-d,.e. (x,y) € S u . 



lu-(x,y) 

The following theorem is proved in JlJ and (|] . 
Theorem 2.1 If there exists a calibration if for u, then u is a Dirichlet minimizer of the Mumford-Shah 



functional (LI) in 



Finally, we present a lemma (proved in Q), which allows to construct a divergence free vector field 
starting from a family of harmonic functions. 

Lemma 2.2 Let U be an open subset of R 2 and I , J be two real intervals. Let u : UxJ — > / be a 
function of class C 1 such that 

• u(-, • ; s) is harmonic for every s S J ; 

• there exists a C 1 function t : Uxl — > J such that u(x, y ; t(x, y ; z)) = z . 
If we define in Uxl the vector field 

4>(x, y, z) := (2Wu(x, y ; t(x, y ; z)), |Vu(x, y ; t(x, y ; z))\ 2 ), 

where \7u(x, y ; t(x, y ; z)) denotes the gradient of u with respect to the variables x,y computed at the 
point (x, y ; t(x, y ; z)) , then (j) is divergence free in Uxl . 

3 Construction of the calibration 

Let {e x ,e v } be the canonical basis in K 2 and for i = 1,2 consider the vectors Tj = (—1/2, (— l) i v / 3/2), 
i/; = ((— l)'v / 3/2, 1/2) , which are tangent and normal to the set Si-i^ (see Fig. 1). As uo(0, 0) < 
iti(0,0) < 1*2(0,0), there exists an open neighbourhood U of (0,0) such that the function u belongs to 
SBV(U), the discontinuity set S u of u on U coincides with {J i<: j(Si.j n U), and the oriented normal 
vector v u to S u is given by 

{vi for (x,y) € S 0> i, 
v 2 for (x,y) € S lj2 , 
e y for (x,y) € ^0,2; 

by the assumptions on , the function u satisfies the Euler conditions for ([hi]) in U . We will construct 
a local calibration ip — (<p xy , <p z ) : UxR — > M 2 xIR for u. 

Applying Schwarz reflection principle with respect to i§q,i and 5*0,2, the function uq can be har- 
monically extended to U \ Si t 2, and analogously u\ and u 2 can be extended to U \ 5o,2 and U \ S*o,i , 
respectively. By the hypothesis on m and by Cauchy-Kowalevski Theorem (see M) the extension of uq 
coincides, up to the sign and to additive constants, with u\ on A\ and with u 2 on A 2 ; analogously, 
the extension of u\ coincides, up to the sign and to an additive constant, with u 2 on A 2 . Since the 
composition of the three reflections with respect to So 1 , $12 ; an, i $0,2 coincides with the reflection with 
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Figure 1: the triple junction. 



respect to the bisecting line of the sector A§ , by the previous remarks we can deduce that uo is either 
symmetric or antisymmetric with respect to the bisecting line of Aq . 

We consider first the case uq symmetric (the antisymmetric case will be studied in Section 7). Then 
also u\ , U2 are symmetric with respect to the bisecting line of A\ , Ai , respectively, and the extensions of 
Uo, Ui, U2 by reflection are well defined and harmonic in the whole set U . 

In order to define the calibration for it, let e > 0, U S (u,_i(0, 0),u,(0, 0)) for i = 1,2, and A > 
be suitable parameters that will be chosen later, and consider the following subsets of UxM.: 

d :— {(x,y,z) £ UxM : Ui(x,y) — e < z < Ui(x,y) + s} fori = 0,1,2, 

Ki := {(x,y,z) EUxR: k + cti(x,y) < z <k + 2\ + (3i(x,y)} fori = 1,2, 
Hi := {(x,y,z) € UxR : k + A/2 < z < k + 3A/2} fori = 1,2, 

where ai and f3i are suitable Lipschitz functions such that ai(0, 0) = /^(O, 0) = 0, which will be defined 
later. If s and A are sufficiently small, then for every i,j the sets Gi , Kj are nonempty and disjoint, 
while for every i the set Hi is compactly contained in Ki , provided U is small enough (see Fig. 2). 

The aim of the definition of the calibration (p in Gi is to provide a divergence free vector field 
satisfying condition (c) and such that 

tp xy (sTi, z) ■ Vi > for < z < Ui^i + e and for ui — e < z < u^, 
tp xy (sTi, z) ■ Vi < for — e < z < and for Ui < z < Ui + e, 

for z = 1 , 2 and s > , and analogously 

if xy (s, 0, z) • e v > for i*o < z < i*o + e and for U2 — e < z < U2, 
(p xv (s, 0, z) ■ e y < for Uq — £ < z < Uo and for U2 < z < «2 + £', 

these properties are crucial in order to obtain (d) and (e) simultaneously. Such a field can be obtained 
by applying the technique shown in Lemma |2.2|, starting from the family of harmonic functions it, + tv. L . 



G 




where we choose as Vi the linear functions defined by 

v {x,y) :=T 2 -{x,y) + e, vi{x,y) := e x -(x,y) + e, v 2 (x,y) := n-(x,y) + s. 
So for every (x, y, z) € Gi , i — 0, 1, 2 , we define the vector ip{x, y, z) as 

z - Ui(x,y), 



n z- Ui (x,y) 

2\/Ui + 2 — VVi, 

Vi{x,y) 



Vi(x,y) 



The role of Ki is to give the exact contribution to the integral in (e). In order to annihilate the 
tangential contribution on S u due to the choice of the field in d , we insert in Ki the region Hi and for 
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every (x, y, z) € Hi, i = 1,2, we define tp(x,y,z) as 

— (Vui_i + Vitj) ,/x 

where ^ is a positive constant which will be suitably chosen later. By the harmonicity of Ui this field 
is divergence free and, as d v ui — on S u for every i , its horizontal component is purely tangential on 
S u . So, it remains to correct only the normal contribution to the integral in (e) due to the field in Gi . 
To realize this purpose on the two segments Si-i,i, i = 1,2, we could require that ai(sTi) — pi(sTi) = 
for every s > (see the definition of Ki ) and define tp(x, y, z) for (x, y, z) € Ki\ Hi as 

^g(Ti-(x,y))vi,(ij , (3.1) 

where g is a function of real variable chosen in such a way that (e) is satisfied for [x, y) € S;_i,i , i.e., 

g(t) := 1 - V* G K, 



as we will see later in ( 3.1G ). Note that the two-dimensional field g {ji-ix, y)) Vi is divergence free, since 
it is with respect to the orthonormal basis {rj, Pi} , hence tp is divergence free in Ki\Hi] moreover, since 
Lp z = /i on Ki, the normal component of ip is continuous across the boundary of Hi , so that tp turns 
out to be divergence free in the sense of distribut ions in the whole set Ki . Actually it is crucial to add a 
component along the direction n to the field in ( ft.l| ) in order to make (d) true, as it will be clear in the 
proof of Step 2 (see Section 5); this component has to be chosen in such a way that it is zero on Si—n 
(so that (e) remains valid on these segments) and that it depends only on • (x, y) (so that the field 



remains divergence free). Therefore we replace in (3.1) the vector g (r, • [x , y)) Vi by 



ct>i(x,y) := {-l) i+1 f {v v {x,y))n+ g{Ti-{x,y))vi, (3.2) 
where / is an even smooth function of real variable such that /(0) = and which will be chosen later 



in a suitable way (see (5.14)). From this definition it follows that 

( f>%(x,y) = - ( f>f(x,-y), $(x,y) = <fi(x,-v), (3.3) 

so that 

fa(x,0) + <h(x,0) =2<f>l(x,0)ey, 
i.e., if we assume that cei(x,0) — (3i(x,0) for every x > 0, the contribution given by the fields (|3.2|) to 



the integral in (e) computed at a point of So, 2 is purely normal, as required in (e), but its modulus is 
in general different from what we need to obtain exactly the normal vector e y . In order to correct it, 
we multiply 4>i by a function Oi which is first defined on Si-ij U So, 2 (more precisely, Oi is taken equal 
to 1 on Si_i,i and to the correcting factor on So, 2); then, we extend it to a neighbourhood of (0,0) by 
assuming Gi constant along the integral curves of cf>i , so that (Jitfii remains divergence free. 

The integral curves of <pi can be represented as the curves {(x, y) £ U : y = ipi(x, s)} , where tpi(x, s) 
is the solution of the problem 

\d x tpi(x,s)(f>f(x, ipi(x,s)) - $(x,ipi(x,s)) = 0, 
[ ipi(s, s) = 0, 
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which is defined in a sufficiently small neighbourhood of (0, 0) . By applying the Implicit Function 
Theorem, it is easy to see that if U is small enough, then there exists a unique smooth function hi 
defined in U such that 

ft<(0,0) = 0, Mx,h l (x,y)) = y. (3.5) 

Note that the curve {(x, y) € U : hi{x, y) = s} coincides with the integral curve {(x, y) € U : y = ipi(x, s)} 
and that (hi(x, y), 0) gives the intersection point of the integral curve passing through (x,y) with the 
x-axis; in other words, the level lines of hi provide a different representation of the integral curves of fa 
in terms of their intersection point with the cc-axis. 
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Figure 3: integral curves of the field fa . 



We state here some properties of hi and ipi for further references. Since ifii(s, s) = 0, we have that 

hi{s,Q)=s (3.6) 



for every s such that (s, 0) G U . By (3.4) and by differentiating the initial condition in (3.4) with respect 
to s , we obtain 

d *MW) = %F77r£ = ^ = L ^r, 9 S ^(0,0) = -9 X ^(0,0) = ^— L-. (3.7) 



€(0,0) if V3 



By differentiating the equation in (3.4) with respect to x and to s, and by using fl3.2| ), it is easy to see 
that 

3tyi(0,0) = Oi(0.0) = 0, (3.8) 
while by differentiating twice with respect to s the initial condition tpi(s, s) — 0, we obtain that 

a s 2 ^(0, 0) = -20,(0, 0) = 0. (3.9) 



By ( |3.7| ) and (|3.8|), the curve {/ij = 0} (which coincides with {y = ipi(x,0)}) is tangent to Ui at 0, 
which may be an inflection point. Moreover, since d x tpi(0, 0) ^ 0, by continuity the function ipi(-,s) is 
strictly monotone in a small neighbourhood of for s sufficiently small; by this fact and by comparing 
the values of the function ipi(-, hi(sTi)) at the points hi(sTi) and srf , it is easy to see that 

hi(sTi) < (3.10) 

for every s > such that sT. t E U , provided U is small enough. Remark that by ([Jl]) and ( 3.1 0| ) it 
follows that the segment So, 2 is all contained in the region {hi > 0} , while Si-i^ in the region {hi < 0} . 
At last, we set 

'1 ithi(x,y)<0, 
<Ti(x,y) := { g(hi(x,y)) 



2</>V(h i (x,y),0) 



if K{x,y) > 0; 



since by definition <^(0,0) = g(0)i/f — g(0)/2, the function Ui is continuous across the curve {hi = 0}. 
Moreover, remark that from (|3.3|) it follows that tp2{x, s) — —ipi(x, s) , h2(x, y) — hi(x, —y) , and then 



For every {x, y, z) € Ki \ Hi 



(X 2 (x,y) = cri(x,-y). 
1,2, we define ip(x,y,z) as 



(3.11) 



In the remaining regions of transition it is convenient to take ip purely vertical. In order to make 
ip divergence free in the whole set UxM., we need the normal component of ip to be continuous across 
the boundary of Gi and Hi . To guarantee this continuity across dGi , we are forced to take as third 
component of ip the function 

|Vu | 2 for z < h + A, 



w(x,y,z) 



e 2 
vl(x,y) 



|Vui| 2 for h + A < z < l 2 + A, 



|Vu 2 | 2 forz>l 2 + X- 



(3.12) 



Finally, wc define the functions a^, [3i in such a way that the normal component of ip turns out to be 
continuous also across the boundary of more precisely, for i = 1,2 we choose cti as the solution of 
the Cauchy problem 

-2 

|Vui_i(x,z/)| 2 , 



1 

-ai(x,y)ct)i(x,y) ■ Vai(x,y) - /i = 5 — 7 \ 



K ai(sTi) = 0, ai(s,0) = for s > 0, 

while /3i as the solution of 

1 e 2 

-a,(x,y)(j)i(x,y) ■ \7(3i(x,y) - /i = ^ r 

A vf(x,y) 

M« T i)=°> A(s,0) = fors>0. 



+ \Vui{x,y)\ 7 
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Since ai is not C 1 near the curve {hi = 0}, we cannot expect a C 1 solution. Nevertheless, if U is 
small enough, then a j , /3j are Lipschitz function defined in U , and the possible discontinuity points of 
Vet,, V/3j concentrate only on the curve {/ij = 0} ; indeed, if U is sufficiently small, the Cauchy problems 



and 



,<3 1 (st 4 ) = (s€R), 
g(hi(x,y)) 



+ IV^-i^.y)! 2 , 



(3.13) 



2A0?(^(a: lV ) > O) 

^(8,0) = o (set), 



<j>i(x,y) ■ V6ii{x,y) - /j, 



\Vui-i(x,y)f 



(3.14) 



admit a unique solution ct^dj € C°°(U), since the lines {sr^ : s G R} and {(s, 0) : s £ R} are 
not characteristic for these equations. Since the curve {hi — 0}, which coincides with the curve {y — 
ipi(x, 0)} , is a characteristic line of both equations ( 3.13 ) and ( ft.14 ) (use ( |3.4| ) and <?(O)/(2A0f (0, 0)) = 1 ), 
the functions dj, d^ assume the same value on the curve {hi = 0} . So, ai can be regarded as the function 
defined by 

&i(x,y) ifhi(x,y)<0, 
&i(x,y) if hi(x,y) > 0, 



and therefore on is C 00 in U \ {hi = 0}, and all derivatives of a, have finite limits on both sides of 
{hi = 0} . The same argument works for /3, . 

The complete definition of the field is therefore the following: for every (x,y, z) £ C/xR, the vector 
(p(x,y,z) = (ip xy 7 <p z )(x,y, z) £ R 2 xR is given by 



vi(x,y) 



Vui 



z - Uj(x,y) 
Vi(x,y) 



(Ti(x,y)<pi(x,y),n 



-y (VUi-l + VUi) ,/x 



mG l (i = 0,1,2), 



in K^H, (t = l,2), 



^ (z = 1,2), 



otherwise. 



By construction conditions (a) and (c) are satisfied. 
Condition (b) is trivial in for all i . 

Since Vitj(0, 0) = for all i (this fact easily follows by the assumptions on the regularity of Ui and 
by the Euler conditions), we have that 



fM" |Vtti(0 '° )|a = 1>0; 
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then, if U is small enough, 



vf(x,y) 



\V Ui (x,y)\ 2 >0 



for every (x, y) G U and for every i = 0, 1,2, and so w is always positive. 

Arguing in a similar way, if we impose that fi > 1/(4A 2 ), condition (b) holds in Ki, provided U is 
sufficiently small. 

By direct computations we find that for every (x, y) € U 



ip xy dz 



e e 
Vvi-i Vvi 

Vi-l Vi 



(fa - a t + X)ai<j)i 



(3.15) 



for i = 1, 2, while 



Note that for i= 1,2 



f U2 e 1 e 2 1 2 

/ ip xy dz = —Wv Vu 2 + T >(& - a t + X)ai(j)i. 

Ju «o «a 



«i_l(sTi) = U»(sTj) = uq(s,0) = -- + e Vs € K, 
V«i-i [x, y) - Vwi (a;, y) = V3vi V(x, y) e [/. 



(3.16) 

(3.17) 
(3.18) 



As hi(sTi) < for every s > by ( 3.10| ), we have that e>i(sTj) = 1 for every s > 0, while by definition 
oti(sTi) = Pi(sTi) = 0. From these facts, (3.15), (3.17), ( ft.l§| ), and the definition of </>j, we obtain 



' ' <p Xy (sT h Z) dz = V3^— H + (-l) i+1 /(0)Ti + g( S )^ = „. 



(3.19) 



where the last equality follows from the definition of g and the fact that /(0) = 0. Analogously, by the 
equalities 



Vo(s,0) = v 2 (s,0) VseR, 
Vv (x, y) - Vu 2 (x, y) = V3e y V(x, y) e J7, 

by the definition of a, and ft, and by (3.3), ( |3.11| ), ( 3.16| ), we have 



(3.20) 
(3.21) 



u 2 {s,0) _ 2 

^( S ,0,z)dz = Vs^— 

u o (s,0) V O (S,0) 



V3 



v o (s,0) 



e y + 2 ( T 1 (s,0)^(s,0)e y 



e y +g{s)e y =e y , 



(3.22) 



where the two last equalities follow from (3.6) and from the definition of <j\ and g . So condition (e) is 
satisfied. 

The proof of condition (d) will be split in the next three sections: in Section 4 we prove that condition 
(d) holds if t\ and t 2 belong to suitable neighbourhoods of tij_i(0,0) and Uj(0, 0), respectively; then, 
in Section 5 we prove condition (d) for t\ and t 2 belonging to suitable neighbourhoods of uo(0, 0) and 
U2(0,0), respectively; finally, in Section 6, by a continuity argument we show that condition (d) is true 
in all other cases. 
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4 Estimates for t\ and ti near U{-\ and U{ 



For (.x, y) G U and ti, t 2 € R, we set 



I(x,y,h,t 2 ) ■= tp xy (x,y,z)dz 



(4.1) 



ti 



and we denote its absolute value by p. In this section, we will show that p(x,y,ti,t 2 ) < 1 in a neigh- 
bourhood of the point (0, 0, Uj_i(0, 0), u,(0, 0)) for i = 1,2, so that the following step will be proved. 

Step 1.- For a suitable choice of the parameter e, there exists S > such that condition (d) holds for 
\t\ — Uj_i(0, 0)| < <5, |t 2 — Wj(0, 0)| < (5 with i = 1,2, provided [/ is small enough. 

Note that p is a continuous function, but its derivatives with respect to x, y may be discontinuous 
at the points (x, y, ti, t 2 ) such that h\(x,y) = or h 2 (x,y) = 0; indeed, the curve {hi = 0} is the 
boundary of the different regions of definition of the functions cr, , on, and /?j , whose derivatives may 
present therefore some discontinuities. Nevertheless, if we set Ni :— {(x,y) e U : hi(x,y) < 0} and 
Pi := {(x,y) G U : hi(x,y) > 0}, the restrictions of <7j, , and to the sets iVj and Pi can be 
extended up to the boundary {hi — 0} as C°° functions; so, along the curve {hi = 0} the traces of the 
derivatives of <7j , on , and are defined. Then, also the traces of the derivatives of p with respect to 
x,y are defined at the points (x,y,ti,t 2 ) with h\{x,y) — or h 2 {x,y) = 0. 



Since we want to study the behaviour of p in a neighbourhood of (0, 0, itj_i(0, 0), u,(0, 0)) , we can 
suppose \t\ — Uj_i(0, 0)| < e and \t 2 — Uj(0, 0)| < £ , so that the possible discontinuities of the derivatives 
of p concentrate only on the curve {/i, = 0} . We study separately the two regions Ni and Pj . 

Consider first the case (x, y) £ Ni, which is the region containing S^-i,, . We will study the derivatives 
of p at the points of the form 



y 




{h,=0| 



Figure 4: the regions Pi and N± . 



qi(s) := (sn,Ui-i(sTi),Ui(sn)), s > 0. 
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We have already shown (condition (e)) that p(qi(s)) = 1 for every s > 0; we want to prove that 

Vp(<7i0)) = Vs>0 (4.2) 

(where now V denotes the gradient with respect to x, y, ti,ta) an d that the Hessian matrix of p with 
respect to v^t\^t<x is negative definite at qi{0). 

Let I Ti and I Ui be the components of the integral in (4.1) along the directions n and v. L , respectively. 
Since by definition 

p(x, y, h,t 2 ) = [{I Ti {x, y, t u t 2 )f + (x, y, hM)?] l/2 , 
the gradient of p is given by 



Vp = -(/ T< v/ Ti +r < vr i ). 
p 



(4.3) 



Note that (3. IE) implies that 

r*fe(s))=0 and = 1 Vs>0, (4.4) 

hence 

Vpfe(s))=Vr*fe(s))- (4.5) 
By the definition of 95 in G,; and by ( 3.1 5| ) we can compute explicitly the expression of I v * at (x, y, t\, t%) ■ 



_ _2{t 1 - Ui-\)d Vi Ui-i + 2(*2 - Ui)d Ui Ui + -{Pi - at + A)<7j 



^ (e 2 - (i! - u^) 2 ) + Ae 2 - (ta - Ml ) 2 ) 



2Vi-! 



2vi 



where 



ti i (x,y) = (-iy +1 f(v i -(x,y)) and ^^y) = g(n-(x,y)). 



(4.6) 



(4.7) 



By differentiating (4.6) with respect to the direction vi we obtain 

d Vi I Vi = 2{d Vl u l _ l ) 2 - 2{d^ Ui ) 2 - 2(ti - Ul _i)a 2 iUl _i + 2(t 2 - u^dl.m 

+ Tpviifii ~ oiijaicj)"* + -{fii - on + \){d Vi ai(/)" i + aiduffl) 

- ^-y - (*i - «-i) 2 ) + ^2 (e 2 - (*a - ^) 2 ) 

H (ti - Ui-i)d Vi Ui-i H (*2 - u l )d v% u l . (4.8) 

By the Euler conditions, 9 Vi ti$_i(sri) = d Vi Ui{sTi) = for every s > 0. Moreover, since |Vttj_i| = |Vu;| 
on [/ (see the remark at the beginning of the proof), in the region iVj the function /3j — ctj coincides 
with the solution £j of the problem 



,6(st;) = (s>0). 



e 2 e 2 



,2 ' 



(4.9) 
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As d Ti ^i(sTi) = and Vi-i(sTi) = Vi(sTi) for every s > (see (3.17)), we have that 

d Vi (Pi - Oi)(sTi) = 9vM*Ti) = °- ( 4 - 10 ) 

By defin ition d Vi 4>^ = and <Ji(x, y) = 1 for every (x, y) € Ni ; using these remarks and the first equality 
in ( 3.17 ) , we can deduce that 



9 Vi r*(q i (a)) = 



(4.11) 



for every s > 0, and the equality holds also for the trace of d Vi l Vi at qi(0) . Since the derivatives of I Vi 
with respect to t\ and t% are given by 



v^3 V3 
d t J Vi = -2d Vi Ui^i (t\ - d t2 I" 1 = 2d Vi Ui (t 2 - Ui), 

Vi-l Vi 

by the Euler conditions it follows that 

d ti n( qi (s))=d t2 I Vi (qi(s))=0. 



(4.12) 



(4.13) 



As I Ui (qi(s)) = 1 for every s > 0, equalities (4.13) imply that d r J Vi (qi(s)) = 0. By this fact, 
( |4.1lD , and ( 4.13[ ), equality (O) is proved. 

Now we need to compute the trace of the Hessian matrix of p with respect to Vi,t\, t 2 at the point 



(4.14) 



<7i(0); using ( |4.4| ) (4.11), (4.13) and (4.2), the Hessian matrix at qi(0) reduces to 

v^ lt2 pfe(o)) = [v^ li2 i Ti ® ^t lt2 i n +Kt l tJ Vi ](n(o)), 



where V„ i t 1 t 2 denotes the gradient with respect to v%.,t\,t% and ® the tensor product. As before, we 
know the explicit expression of I Ti : 

I Ti = -2(ti - Ui-i)d n Ui-i + 2(t 2 - Ui)d Tz Ui + -(Pi - a t + A)er 4 0,[ l 

A 

- ^—(e 2 ~ (h - u^f ) + ^-(e 2 - (t 2 utf), (4.15) 
2vi-i 2vi 



hence, using the Euler conditions, (4.10), and the fact that Oi = 1 in Ni, it results that 

^r*te(o)) = l dviVi _ 1 (o ! o)-l dvi v i (^o) + 9^(0,0) = ^, 



where the last equality follows by ( 3. 18| ) and by the equality 

d Ui <j>?(o) = (-i) i+1 f(o) = o. 



(4.16) 



(4.17) 



By differentiating (4.8) and by using the Euler conditions, ( 4.1 0[ ) , the constancy of ct; in Ni, and the 
fact that dl.cj)^ = 0, we have 

dlJ Vi (qm) = ^(0,0)^(^-^(0,0) + ^d Vi Vi^(0,0) - ^9^(0,0) = (4.18) 
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where the last equality follows from 

1 



2\/3 



C (0,0)9^ (A - o<)(0,0) = -— , (4.19) 

A S 

which can be obtained by differentiating ( f4.9| ). Using (4.14), (4.1€), and ( |4.18j ), we obtain that 

= [^fe(O))] 2 + ^r*fe(0)) = | - ^ < 0, (4.20) 

provided e is sufficiently small. Since dt 1 I T '(qi(0)) — (this can be easily proved using the fact that 
Vui_i(0, 0) = Vu 4 (0, 0) = 0), by (glj) we have that 



By differentiating ( 4.12 ) and by using the Euler conditions, it turns out that 



3Ur'(?<(0)) = -2# 4 u<_i(0,0), 



so that 



det 



dl ltl P dip 
Arguing in a similar way, one can find that 



(*(0)) = o_2 



3 

2? 



1 - 



dll Ui (qi(0)) 
V3 ' 



e 



e\ -4«^_r(0,0)r. 



dl it2 p(u(°)) = 2dlM°> °). a? a P(»(0)) = 9 t 2 i42 p( % (0)) = 0, 



so that 



detV^ it2(O (g z (0)) = - 



3%/3 
2^ 



4^3 



(^Ui-iCO.O)) 2 + 0)) a 



Since for e sufficiently small it results that 



det 



dZ.P du itl P 
9i itl P dip 



( 9i (0)) >0, detV* tlta p(eft(0))<0, 



(4.21) 



then, by ( 4.20 ) and ( 4.21 ) the Hessian matrix of p at qi(0) is negative definite. 

At this point we have all the ingredients we need in order to compare the value of p on Si—i,, with 
its value at a point (x,y,t%,t2) for (x,y) G AT, and |ti — Uj_i(0, 0)| < e, |t2 — ih(0, 0)| < e. 

Remark that since the curve {hi = 0} may have an inflection point at the origin, the set Ni might 
be not convex. If the segment joining (x, y) with its orthogonal projection on Si-ij (which is a point 
of the form sTi with s > ) is all contained in Ni , then we can consider the restriction of p to the 
segment joining (x,y,ti,t2) with qi(s) and write its Taylor expansion of second order centred at qi(s) . 
By (4.2) and the fact that the Hessian matrix of p is negative definite at qi(0) (and then, by continuity 
in a small neighbourhood), we have that there exist S, C > such that, if U is small enough and 
|ti -Ui_i(0, 0)| < 5, 1*2-^(0,0)1 < S, then 



p(x,y,h,t 2 ) <1 - Civ.-ix^y)) 2 - C{h - u^isTi)) 2 - C{t 2 ~ Ul (sTi)) 2 < 1. 



1G 



In the general case, since the curve {y = i/)i(x,0)} is C 2 with null second derivative at 0, one can find 
s > , a G K such that the segment joining {x, y) with STi + ai>i is all contained in Ni and the ratio 
\(x,y) — sTi — avi\/ci 2 is infinitesimal as a — > 0. Since s > 0, the segment joining STi + avi with its 
projection STi on Si-% t i is all contained in iVj, so that we can apply to this point the estimate above; if 
we call L the L 00 -norm of the gradient of p, we obtain that 

p(x,y,ti,t 2 ) < p(sTi + avi,ti,h) +L\(x,y)- 8Ti- aui\ 



< 1 — a [ C — L 



(x,y)- sn - avi\ 



- C{h - U^STi)) 2 ~ C(t 2 - U^STi)) 2 



which is less than 1, provided U is small enough. So we have proved that, if e is sufficiently small, then 
there exists S > such that 



p(x,y,t 1 ,t 2 ) <l for (x,y) <G Ni, \t x - Ui_i(0, 0)| < 6, \t 2 - Ui(0, 0)| < 5, 



(4.22) 



provided U is sufficiently small. 

Suppose now (x,y) € P i: \ti — u,_i(0,0)| < e, \tz — Wj(0,0)| < e. In order to show that p < 1 
also in this case, we will compute the traces of the gradient and of the Hessian matrix of p at the point 
qi(0) . The main difference with respect to the previous case is that in the region Pi the function /3j — oti 
coincides with the solution r/i of the problem 



vt_ x {x,y) vf{x,y) 



(4.23) 



.77*0,0) =0 (s>0), 
while the function Oi is defined as 

g(hi(x,y)) 



&i(x,y) 



2$(h i (x,y),0) 



V(x,y)EP l . 



By (4.4) and (Hj) it follows that 



Vp(«(0)) = Vf i (? j (0)). 



(4.24) 



(4.25) 



By (4^6) we obtain the following expression for the gradient of I Vi with respect to Tj, Vi computed at the 
point qi(0): 



V T4! , i r i (gi(0)) = fl (0)V«7,(0,0) + Vf(0,0) + ^-n, 



(4.26) 



where we have used the Euler conditions, the fact that V(/3, — aj)(0, 0) = by (4.23), and that 

V«i_i (x, y) + Wvi (x, y) = -n V(x, y) G U. 



It follows immediately by (4.7) that 



V(j)i*(x,y) = g'(n-(x,y))i 



(4.27) 
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and by the definition of g that 



«g(i,o) 



«o 2 (^0) 



for all teR. By differentiating ( 4.24 ), we obtain that 



Vai(x,y) = -p(hi(x,y))Vhi(x,y), 



(4.28) 



(4.29) 



where we have set 



J>(*) 



g'(t) g(t) 

^(t,0) [#(i,0)] 



r&0f(i,O). 



To compute the gradient of /ij it is enough to differentiate the second equality in (3.5): this provides 

d x ipi(x, hi) + d a ipi(x, hi)d x hi = 0, d s ipi(x, hi)d y h t = 1; (4.30) 



by (3.7) we have that 



Since 



V/ij(0,0) 



-2r,: 



(4.31) 



fl^Or, y) = {-ly+^-fiyi ■ (x, y)) - \g'{ Ti ■ (z, y)), 



we find that p(0) = 3g'(0)/g(0) , and substituting in (4.29), we have that 

V.,(0,0) = -3^|„. 



(4.32) 



Since the partial derivatives of I v * with respect to t\ and £2 are still given by (4.12), they are equal to 
at the point qi(0) , as in the previous case. Then, by ( ggg ), ( pL26| ), ( ggg ), (|4.32j), and ( ggg ), we deduce 
that 



Vp(ft(0)) 



3^3 



TV, 0,0 



(4.33) 



To conclude the study of p in this region, we write the Hessian matrix of p with respect to Vi,t\,t%, which 
still satisfies (4.14). Differentiating (4.15) and using the Euler conditions, the fact that V(/3, — ccj)(0, 0) = 
0, 0^(0,0) = and (4.17), we obtain that (4.16) still holds. Differentiating (4.8) and computing the 
result at 9i(0), we have that 

^te(O)) = ^(0)^(ft - Oi)(0,0) +$(0)3^,(0,0) + l(ft /i « i _ 1 (0,0) - 9^(0,0)), (4.34) 

where we have used in particular that 5^(7,(0, 0) = by ( |4.32| ) and that d^cj)"* = 0. In order to compute 
the second derivative of pi — on with respect to the direction Vi , we differentiate ( [4.25 ) with respect to 
x and with respect to y; using the fact that 9 x (ft — «i)(s,0) = for every s > 0, we obtain 



^(ft - Oi)(0, 0) = 0, dUfc - 00(0,0) = -(-l) i+1 — 



^(ft"«i)(0,0) 



2V3 A 
£ 5(0) 



+ V3(-ir +i C(A- ai )(o,o) = 



.9(0)' 
4^3 A 



£ 5(0)' 



(4.35) 
(4.36) 
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By the relation d*. = \d 2 x + ^(-l) l dl y + \d 2 y , it follows that 



c£ (A-a,)(0,0) 
Since 9^(0,0) = by ( |4.3lj) , from (|^|) wc obtain that 



5(0) 



5'(0) 



9(0) 



3fl'(0) 



0^(0,0)) <fti(0,0) = -^c^(0) 



>r(o,o) [^(o,o)] 2 

By differentiating twice with respect to the direction vi the second equality in ([T^), we obtain that 

{v^fdl^ix, hi) + 2is?d 2 xs Mx, hi)d Wt h t + d 2 i>i(x, hi)(d Ui hi) 2 + d s Mx, hi)d 2 .hi = 0; 



(4.37) 



since 3„ i,(0 0) = by ( |4.31[ ) and 9^(0,0) = by (|J|), we can c onclu de th at dg > t (0,0) = and 
then, by (4.37) also the limit ol d 2 .Oi at (0,0) is equal to 0. Taking ( [3.17 ) and (4.34) into account, we 
can conclude that 

2e 



i.e., (4.18) is still satisfied. Since it is easy to see that also the other second derivatives of p remain 
unchanged, we can conclude that the Hessian matrix of p with respect to Vi,t\, £2 is negative definite at 
<7i(0). _ 

If the segment joining (x, J/, tijta) with %(0) is all contained in Pi , then we consider the Taylor ex- 
pansion of second order centred at qi{0) of the funct ion p restricted to this segment; since the component 
of (x,y) along is less or equal than 0, by ( 4.33 ) and by the fact that the Hessian matrix of p with 
respect to v%,t\,t2 is negative definite, we have that there exists S > such that p(x,y,ti,t2) < 1 for 
\t% — u-i— 1 (0, 0) I < S, \t2 — Uj(0,0)| < S, provided U is small enough. In the general case, we can find 
s < 0, a G K such that the segments joining (x,y) with sr; + Wi, and sTi + avi with (0,0) are all 
contained in Pi , and \{x, y) — STi — aVi\/a 2 is infinitesimal as a — > 0. Arguing as for the region iV, , this 
is enough to obtain the same conclusion. So we have proved that, if e is small enough, there exists S > 
such that 



p(x,y,tx,t 2 ) < 1 
provided U is sufficiently small. 



for fay) e Pi, |*i - itt-i(0,0)| < 5, |*2-u<(0,0)| < 6, 



(4.38) 



By (4.22) and (13§) Step 1 is proved. 



□ 



5 Estimates for t\ and near uq and 112 

This section is devoted to the proof of the following step. 

Step 2.- For a suitable choice of the function / (see ( |3.2| )), there exists 5 > such that condition (d) 
holds for \ti — uq(0, 0)| < 5, | *2 — "2(0, 0)| < S, provided U is small enough. 

In order to prove the step, we want to show that the function p, introduced at the beginning of Section 4, 
is less or equal than 1 in a neighbourhood of the point (0, 0, «o(0, 0), ^2(0, 0)) . We can assume that 
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\ti — uo(0, 0)| < e, \t 2 —1*2(0, 0)| < e. Since now the derivatives of p may be discontinuous on the curves 
{hi = 0} and {h 2 — 0}, we have to consider separately four different cases, one for (x,y) belonging to 
each one of the regions N x n N 2 , N x n P 2 , N 2 C\P 1 , and P 1 C\P 2 . 

Let P and I v be the components of the integral in (4.1) with respect to e x and e v , that are the 
tangent and the normal direction, respectively, to the third part of the discontinuity set So j2 . 

Consider first the case (x, y) € Pi fl P 2 , which is the region containing So,2 ; as before, we will study 
the derivatives of p at the points of the form 



q (x) := (x,0,u (x,0),u 2 (x,0)), 



x>0. 



Condition ( 3.22 ) implies that p(qo(x)) = 1 for every x > 0; we want to prove that 

\7p(q {x)) = Mx > 



(5.1) 



and that the Hessian matrix of p with respect to y,t±,t 2 is negative definite at Qi(0) . By the definition 
of p , it follows that 

vp = -(rwr + I V VP). 

P 

Since I x (qo(x)) = and P(qo(x)) = 1 for every x > 0, we have that 

Vp(q (x))=VP(q (x)). 

By ( |3.16 ) and by the definition of ip in G,: we can write the explicit expression of P at the point 

{x,y,h,t 2 ): 



1 2 

7 s = -2(<i - u Q )d y UQ + 2(t 2 - u 2 )d y u 2 + - J^(/3 4 - Q, + X)ai<f^ 



2v 



V3 
2v 2 



(s 2 (t 2 



u 2 ) 



(5.2) 



and by differentiating with respect to y , we obtain 

d y P = 2(d y u f - 2{d y u 2 f - 2(ti - u )d> + 2(t 2 - u 2 )9> 2 



1 z Q 

3 



2=1 



«2) 



(*i - uo)d y ua + ~^(h - u 2 )d y u 2 . (5.3) 



Since in the region P x n P 2 the functions /3j — oti coincide with the solutions of the problems (4.23), it 
results that d y (/3i — (Xi)(x, 0) = for i = 1,2. Moreover, differentiating (3.11) and the second equality 
in (3.3) with respect to y, we have that 



d v (J 2 {x,y) = -d y ax{x,-y), d y <f%(x,y) = -d y <j>\{x,-y), 



(5.4) 



and then, using again ( J3.3| ) and (3.11), 

4>\{x, 0)d y a 1 (x, 0) = -<f%(x, 0)d y a 2 (x, 0), a^x, Q)d y 4>\{x, 0) = -a 2 (x, 0)dy<f%(x, 0). 
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By the Euler conditions, d y uo(x, 0) = d y ui{x, 0) = for every x > 0; using all these remarks and ( 3.2CQ , 
we deduce that d y I y (qo(x)) = for every x > and the equality holds also for the trace of d y I v at 
<7o(0) ■ Since we have that 

fo fo 

d tl P = -2d y u - — (ti-uo), d t2 P = 1d y u% - — (fa - U2), (5.5) 

by the Euler conditions it follows that d tl I v (qo(x)) — d t2 I y (qo{x)) = 0. As T y (q (xj) = 1 for every 
x > 0, this implies that d x I v (qo(x)) = 0. Thus we have obtained equality (5.1). 
By (5.1) and (3.22) the Hessian matrix of p computed at <7q(0) reduces to 

V^ li2 p(<Zo(0)) = [VyutJ* ®V vtl t 2 I x + V ytlt2 I y ](q O (0)). (5.6) 
As before, we know that 

1 2 

l x = -2(t a - u )d x u + 2(t 2 - u 2 )d x u 2 + - 22(Pi -a l + Xja^f 

" " (il " " o)2) + 2^ (£2 " (<2 " U2)2) ' 

hence, by the Euler condition, the fact that d y {(3i — a,)(0, 0) = for i = 1,2, and ( 3.21 ), it results that 

V*(?o(0)) - ^+X)$,M?)(Q,0) - ^+2^(0,0) +2^(0,0)^(0,0), 

i=l 

first equalities in (3J3) and in (5.4), and the relation d v 4> x (x, y) = d y <pf(x, —y) . 



where we ha ve also used the 
From ( 4.32 ) we obtain that 



^(0,0) = — m 



Then, using the definition of (f>f and (4.28), we can conclude that 



9^(0,0) = ^-3.g'(0)= 2^3. 



(5.7) 



By differentiating (5.3) with respect to y and by using the Euler condition and the fact that d y (/3i 
cti) (0,0) = for i = 1,2, we obtain 



d 2 y I V (qo(0)) 



Equality ( 4.36| ) implies that 



j - chW + <£M?)](0, 0) + 

2—1 



l^[^(ft-a i )<r^](0,0) = ^. 

7.--1 



(5.8) 
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In order to write explicitly dpi at (0,0), we differentiate the y-component in (4.29) with respect to y 
and we pass to the limit, taking into account that d y hi(0) = (— l) t+1 y/3 by (4.31): 

flgffi(0 > 0) = |p'(0) + ip(0)^ft i (Q). 
By differentiating with respect to y the second equality in ( 4.30| ), we obtain that 



dlhi (0,0) 



(8 h rn n ^ ^i(0,0) 
' (9 ^ l(0 ' 0)) a^o) 



o. 



where the last equality follows by (^9). Since 

,.9"(0) 



p'(0) 



5(0) 9 2 (0) 



[g'(0)] 2 ^jf(0,0) 
3(0) 



while 



«(0,0) 



3^3 ,„ 



/"(0) + - 3 "(0), ##(0,0) 



~r(Q) 



fl"(0), 



and <?"(0) = — v / 3/(2e) , we can write that 



i=i 



y^OS i -a i + A)flg(ff i ^)(0,0) = (20^(7! +4^,0-!^ + 2fig#) (0,0) 

= 2V3/"(0) 



2V3/"(0) + 3 ff "(0) 
3^ 



2e 



Substituting (5.8) and ( 5.11 ) in the expression of dyl y , we find that 



dpy(q o (0)) = 2V3/"(0) 
From ( |5.6| ), (|5.7|), and ( 5.12| ), we finally obtain that 



4^3 



9>( 90 (0)) = [d y I x (q (0))} 2 + d 2 y iy(q (0)) = 12 



4V3 



2 v / 3/"(0). 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



As in the previous step, we can compute explicitly the other elements of the Hessian matrix of p and we 
find that 



dot 



d yP d y tl P 
dytiP d^p 



M0)) 



detV 2 tit2 p(g (0)) 



6^3 ,„ 



r(o) 



-/"(0)-^-^-4(9> (0,0)) 2 , 
- ^ + l^p> (0,0)) 2 + (9> 2 (0,0)) 2 ] 



If we impose the following condition on the second derivative of / at : 



/"(0) < -2V3- 



2s 



(9 2 u (0,0)) 2 + (9 2 u 2 (0,0)) 2 



(5.14) 
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then the Hessian matrix of p is negative definite at qo(0)- 

To conclude, we restrict p to the segment joining (x, y,t\,t 2 ) with qo(x) and we write its Taylor 
expansion of second order centred at q (x) ; using ( |5.1[ ) and choosing / satisfying ( [3.14 ) (so that the 
Hessian matrix of p is negative definite at <?o(0) , and then by continuity in a small neighbourhood), we 
obtain that there exists 6 > such that 



p(x,y,t u t 2 ) < 1 for [x,y) £P 1 nP 2 , \h - u„(0, 0)| < 6, \t 2 - u 2 (0, 0)| < S, 



(5.15) 



provided U is sufficiently small. 

Let us consider the set N\ PI N 2 ■ in this region a\ = a 2 = 1 , while the functions /3j — coincide 



with the solutions of the problems (|4.9|). By (3.22) the gradient of p at the point qo(0) is given by 

Vp(«,(0)) - VP(go(0)). (5.16) 



By ( |5.2D we derive the explicit expression for the gradient of I v with respect to x,y; using the Euler 
condition, the fact that V(/3, — «i)(0,0) = 0, the constancy of Oi and the equality 



Vv (x,y) +Wv 2 (x,y) = 



V(x,y)EU, 



(5.17) 



we obtain that 



V3 
2 



3V3 



Since the partial derivatives of I y with respect to t\ and t2 are still given by (5.5), they are equal to 
at qo(0) , as in the previous case. Therefore, we have that 



VpMo)) 



(5.18) 



If {x,y) 7^ (0,0) belongs to Ni D N 2 and the segment joining (x,y) with (0,0) is all contained in 
NiSlN 2 , then by the Mean Value Theorem, (5.18) and the fact that x is strictly negative, we can conclude 
that there exists 5 > such that 



p(x,y,ti,t 2 ) < 1 for |ti — ua (0, 0) | < 5, \t 2 - u 2 (0, 0)| < 8, 



(5.19) 



provided U is sufficiently small. If the segment joining (x, y) with (0, 0) is not contained in Nir)N 2 , then 
we can find a regular curve connecting (x, y) and (0, 0) , along which we can repeat the same estimate 
as above. 

At last consider the set N 2 n P\ , since the case N\ n P 2 is completely analogous. In this region, <j\ 
is defined by (4.24), while a 2 is identically equal to 1; the function (3i — a± coincides with the solution 



of the problem (4.23) for i = 1, while f3 2 — a 2 with the one of (4.9) for i = 2. Equality (5.16) still holds, 
as well as the fact that V(/3j — o*i)(0, 0) = (0,0) for all i; since Vci is given by the formula (4.32) and 
Vct 2 = 0, by (HH), (|J), and ( |t| ) we have that 



V xy I y (q (0)) 



2 

E 



V^(0,0) 



6^(0,0^(0,0) + ^e 2 
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hence 



VpMo)) 



3\/3 



r 2 ,0,0 



Since the gradient of p vanishes along the direction [yi, 0, 0) , we need to compute the Hessian matrix of 
p with respect to V2,ti,t 2 at the point qo(0); from the equality Vi /2tlt2 / i '((7o(0)) = 0, we have that 

V' 2tlt2 pM0)) = [Vv 2tl t 2 I x ® V V2tlt2 P + Vl 2tit jv](q (0)). (5.20) 

Using the fact that Vu o (0, 0) = Vu 2 (0, 0) = and V(/3< - c*i)(0, 0) = 0, we obtain 

2 i 

d V2 I x (q (0)) = £ft s #(0,0) + 3^(0,0)^(0,0) + -d V2 (v - v 2 ) 
i=l 

= a u ^(o,o)-^'(o) + ^ = >/3, 

where the second equality follows from ( 4.32 ) and from the fact that d V2 §\ + du 2 ( t ) 2 = dy^i & t (0,0). 
If we differentiate (5.2) twice with respect to the direction vi and we compute the result at the point 
<7o(0) , we obtain 



d 2 V2 F (0, 0) = ( i Yl 9 i (A - «<) ff *^ + E M + + 2a,oift*#) (0, 0) + 

\ t=l 8=1 / 

From (4.35) and (4.36), and from ( 4.19[ ) it follows respectively that 

4%/3 A 



3V3 
IT' 



C(A-ai)(0,0) = 



9* 2 (/3 2 -a 2 )(0,0) = 



2^3 A 



Since by ( 4.29 ) we have that d^ 2 ai(x,y) = hp(hi(x, y))d V2 hi(x, y) , then 

d 2 „M0,0) = ip'(0)(^ 2 /n(0,0)) 2 + ip(Q)c£>(0,0). 
Some easy computations show that d^ 2 h%(0,0) = 0; using ( 4.31 ) it results that 



(5.21) 



(5.22) 



^ 2 <7 1 (0,0) = |p / (0) 



9[5'(0)] 2 , ^/"(O) 



2 g*(0) 



3(0) 



(5.23) 



where the last equality follows by (5.9) and by the first equality in (5.1C). At last, by using (3.3) and 
that 

2 

JX#(0,Q) = 1^^(0,0) + ^#(0,0) = -|V3/"(0) + ^"(0), 



(5.10), we obtain that 

2 



(5.24) 



i=i 



and by substituting ( |5.22| ), ( 5.23 ), and ( 5.24 ) in ( 5.21 ), we deduce that 



d 2 V2 P(q (0)) = V/"(°) + Z T- 
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hence 



^(9o(0))=3+^ + ^/"(0) 



By differentiating ( p.5| ) with respect to v 2 and by (5.2C), we obtain 

3Up(«>(0)) = -2a, 2 9 yM o(0,0) = -9>o(0,0), dl 2t2 p(qo(0)) = 2d V2 d y u 2 (0,0) = d 2 y u 2 (0,0). 
At this point, it is easy to see that, if / satisfies the condition 



/"(0) < -2V3- 



6 



p> (o,o)) 2 + (a> 2 (o,o)) 2 



(5.25) 



then the Hessian matrix of p with respect to ^2,*i,*2 is negative definite at the point go(0). Arguing 
as for the region Pi in the previous section, it can be proved that, if / satisfies ( 5.25 ), then there exists 
S > such that 



p(x,y,t u t 2 )<l foT(x,y)eN 2 nPi, |* x - u (0, 0)| < 5, |* 2 - u 2 (0, 0)| < 6, 



(5.26) 



provided U is sufficiently small. 

Since condition (5.14) implies (5.25), if we require that ( |5l^) holds, then by ( jsl^ , |yj ), and |[2§), 
we can conclude that Step 2 is true. □ 



6 Proof of condition (d) 

In this section we complete the proof of condition (d). To this aim it is enough to check condition (d) in 
the three cases studied in the following step, as it will be clear at the end of the section. 



Step 3.- If e is sufficiently small, S S (0, e) , and U is sufficiently small, condition (d) is true for t\ < t 2 
whenever one of the following three conditions is satisfied: 

1) |*i-tt o (0,Q)| >S and |*i — ui(0,0)| > 5; 

2) |t 2 - wi(0, 0)| >S and |t 2 -« a (0,0)| > S; 

3) 1*1-1*0(0,0)1 > § and \t 2 - u 2 {0, 0)| > 5. 



Let us fix 5 £ (0, e) and set 

Mi(x,y) := max{\I(x,y,U,t 2 )\ : u (x,y) - e <h <t 2 < u 2 (x,y) + s, 

1*1-^0(0,0)1 > 5, 1*1-^(0,0)1 >S}. 

It is easy to see that the function M\ is continuous. Let us prove that Mi (0,0) < 1. For simplicity of 
notation, from now on we will denote 7(0, 0,*i,* 2 ) simply by I{t\,t 2 ) and Ui(0, 0) by 

Let *i,*2 be such that uq — s < t± < t 2 < u 2 + e with |*i — uq\ > 5 and |*i — u\\ > 6. Suppose 
furthermore that |*i — ui| < e; then, we can write 

/(*i,* 2 ) = J(*i,«i) +I(u x ,u 2 ) +I[u 2 ,t 2 ), 

I{u 2 ,t 2 ) = I(u 2l t 2 V (u 2 - e)) + I(u 2 - e,t 2 A (u 2 - e)). 
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Therefore, we have 

I{h,t 2 ) = + I(u 1 ,u 2 ) + I{u 2 ,t 2 V (u 2 - e)) - I(t 2 A (u 2 - e),u 2 - e). (6.1) 
From the definition of ip in Gi,G 2 it follows that 

I{si,ui) = - ui) 2 e x for |s x - u x | < e, 

7(u 2 , s 2 ) = -(s 2 - u 2 ) 2 ti for |s 2 - u 2 | < e; (6.2) 

e 

using condition (e), we have that 

S 2 

I(h, ui) + Z(ui, u 2 ) + /(u 2) t 2 V (u 2 - e)) G ^ 2 e x + (6.3) 

e 

where i?i is the parallelogram spanned by the vectors et\ and — [e — ^-^ e 1 . Let C be the intersection 

of the half-plane {(x, y) £ R 2 : u 2 -{x, y) > 1 — v^e} with the open ball centred at with radius 1; some 
elementary geometric considerations show that 

6 2 

v 2 e x + R X C C. (6.4) 

£ 

If Tj is the segment joining with g(Q)vi, then from the definition of ip in Ki, it follows that 

I(ui-i+e,v*-e) = g(0)vi, (6.5) 

and 

I(s u s 2 )£T i (6.6) 



for + e < Si < s 2 < m — e, i = 1,2. Let D := — T 2 ; from (|6.l[), fl6.3j), (|6.4|), and (|6.6|), we deduce 
that 

I(h,t 2 ) £C + D; 

since g(0) = 1 — %/3e, the set C + D is contained in the open ball centred at with radius 1. This 
concludes the proof when \b\ — u\\ < e. 

If \t 2 — i*i | < e, we consider the decomposition 

I(h,h) = I{h,u ) +I(u ,ui) +I(ui,t 2 ), 

I(t 1 ,u ) = I(ti A (u + e),u ) + I(t 1 V (uq + e),uo + e), 

and the proof is completely analogous. 

When |ti — ui\ > e and |i 2 — u%\ > £ , we can write 

I(ti,te) = I(ti,u ) +I(u ,u 2 ) +I(u 2 ,t 2 ), 

I(h,u ) = 7(ti A (tto + e), uq) + Z(ti V (wq + e), u a + e), 

I(u 2 ,t 2 ) = 7(u 2 , i 2 V (u 2 — e)) + 7(u 2 — e, i 2 A (u 2 — e)); 
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therefore, we have 

I{t\,t 2 ) = I(ti A (ito + e), uo) + ^("o, U2) + I(u 2 ,t 2 V (u 2 - e)) 

+ 7"(ti V (n + e), t 2 A (u 2 - e)) - I(u +e,u 2 - s). (6.7) 

Since from the definition of if in Go it follows that 

7(so,uo) = --(so - u ) 2 t 2 for \s - u \ < e, (6.8) 

£ 



using condition (e) and ( |6.2| ), we have that 

I(ti A (it + e), u ) + 7(-u , u 2 ) + I{u 2 , t 2 V (it 2 - e)) e e y - —r 2 + R 2 , (6.9) 

where R 2 is the parallelogram spanned by the vectors er\ and — i^e — r 2 . Let E 1 be the parallelogram 

having as consecutive sides T\ and T 2 , and let F be the set E — g(0)e v ; as /(ui — e,iti + e) = 0, from 
(S.5) it follows that 

I(u +£,u 2 -e)= g(0)e v = (1 - VSe)e v , (6.10) 

and from ( |6.6| ). 

I(si,s 2 )eE (6.11) 

for every u + e < Si < s 2 < u 2 — e, with |si — u%\ > e and |s 2 — u%\ > e. From ( |6.7[ ), (6.9), d6.1Gj ), 
(6.11), we obtain that 

i"(*i,*a) e e y t 2 + R 2 + F. 

£ 

The set e y — -t 2 + R 2 + F is a polygon, since it is the sum of two polygons, and it is possible to prove 
that, if e < y3, its vertices are all contained in the open ball with centre and radius 1. Then, under 
this condition, the whole set e y — ^-t 2 + R 2 + F is contained in this ball; this concludes the proof of the 
inequality Mi (0,0) < 1. 

By continuity, choosing U small enough, we obtain that M\(x, y) < 1 for every (x,y) € U, which 
proves 1). 

To prove 2) and 3), we define analogously 

M 2 (x,y) := max.{\I(x,y,ti,t 2 )\ : u (x,y) - £ <h <t 2 < u 2 (x,y) +s, 

|*a-«i(0,0)| > S, \t 2 -u 2 {0,0)\>5}, 



M 3 (x,y) := max{\I(x,y,ti,t 2 )\ : u (x,y) - e < ti < t 2 < u 2 (x,y) + e, 

|*i — tio(0,0)| >S, 1*2-^(0, 0)1 >5}. 

It is easy to see that the functions M 2 and M3 are continuous and, arguing as in the case of M\ , we can 
prove that M 2 (0, 0) < 1 and M 3 (0, 0) < 1 , which yield 2) and 3) by continuity. □ 
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Conclusion.- As in Step 3, we simply write Ui instead of Uj(0,0). Let us show that, if / satisfies 



( [5.14 ), and e and U are sufficiently small, then condition (d) is true for u a (x, y)—£ < t\ < t 2 < u 2 (x, y)+e 
and in fact for every t\.t 2 G R, since tp xy (x,y, z) = for z < u (x,y) — e and for z > u 2 (x,y) + e. 

We start by considering the case \ti ~ u \ < S . If |ia — ui\ < S, the conclusion follows from Step 1. If 
1*2 — Mi| > S , the result is a consequence of Step 2 when |t 2 — ^2 1 < 6, and of Step 3.2) in the other case. 

We consider now the case \ti — uq\ > 6. If |<i — Wj.| > S, the conclusion follows from Step 3.1). If 
\ti — Ui\ < 5 , the result is a consequence of Step 1 when \t 2 — u 2 \ < S , and of Step 3.3) in the other case. 



This concludes the proof of condition (d) and then, of Theorem 1.1 in the case uq symmetric. □ 



7 The antisymmetric case 

In this section we show how the construction of the calibration for Ui symmetric can be adapted to the 
antisymmetric case. 

If the function uq is antisymmetric with respect to the bisecting line of Aq , then the reflection of uq 
with respect to the S^i and to 50,2 provides an extension of uq , which is harmonic only on £1 \ Si^ 2 and 
which is multi-valued on Si j2 , since the traces of the tangential derivatives of u$ on S\^ 2 have different 
signs. Since ui 7 u 2 coincide, up to the sign and to additive constants, with the reflections of uq with 
respect to So,i and 6*0,2, respectively, they are antisymmetric with respect to the bisecting line of A± 
and A 2 , respectively, and then, their extensions by reflection are harmonic only on f2\ ^0,2 and fi\ Sq,i , 
respectively. 

The calibration ip can be defined as before, just replacing the sets Go,Gi, G 2 with 

G = {(x,y,z) £ {U\S h2 )xR : u (x,y) 
Gi = {(x,y,z) e (C/\5 , 2 )xK : u x {x,y) 
G 2 = {(x,y,z) E (U\S 0A )xR: u 2 (x,y) 

and the sets Hi , H 2 with 

Hi = {{x,y,z)e(U\(S h2 USo, 2 ))xR:li + X/2<z<l 1 + 3\/2}, 
H 2 = {(x,y,z)e(U\{So.ii)S , 2 ))^-l2 + X/2<z<l 2 + 3X/2}. 

Since uq is harmonic in ft \ Sj.,2, the field ip is divergence free in Go by Lemma ^.2[ Moreover, the 
normal component of tp is continuous across the boundary of Go since d„ 2 uo = d V2 v$ = on Si. 2 . The 
same argument works for the sets Gi,G 2 . By the harmonicity of uo and ui, the field is divergence free 
in Hi and the normal component of tp is continuous across the boundary of H± since d U2 UQ — on Si t2 
and dyUi — on 5o,2 ■ Therefore, condition (a) is still satisfied in the sense of distributions on UxR. 
It is easy to see that conditions (b), (c), and (e) are satisfied. 

The proof of Step 1, Step 2, and Step 3 can be easily adapted; indeed, even if now the function 
\I(x, y, tt, f2)| may present some discontinuities when (x,y) £ Sjj, we can write U as the union of 
finitely many Lipschitz open subsets Ui such that |/| is C 2 (UixM. 2 ) and study the behaviour of |7| 
separately in each Ui . So, it results that also condition (d) is true. □ 



— e < z < uq(x, y) + e}, 
-e < z < ui(x,y) +e}, 

- e < z < u 2 (x, y) + e\, 
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